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ABSTRACT 

We study how the first galaxies were assembled under feedback from the accretion onto a central 
black hole (BH) that is left behind by the first generation of metal-free stars through self-consistent, 
cosmological simulations. X-ray radiation from the accretion of gas onto BH remnants of Population HI 
(Pop III) stars, or from high-mass X-ray binaries (HMXBs), again involving Pop HI stars, influences 
the mode of second generation star formation. We track the evolution of the black hole accretion 
rate and the associated X-ray feedback starting with the death of the Pop HI progenitor star inside a 
minihalo and following the subsequent evolution of the black hole as the minihalo grows to become an 
atomically cooling galaxy. We find that X-ray photoionization heating from a stellar-mass BH is able 
to quench further star formation in the host halo at all times before the halo enters the atomic cooling 
phase. X-ray radiation from a HMXB, assuming a luminosity close to the Eddington value, exerts an 
even stronger, and more diverse, feedback on star formation. It photoheats the gas inside the host 
halo, but also promotes the formation of molecular hydrogen and cooling of gas in the intergalactic 
medium and in nearby minihalos, leading to a net increase in the number of stars formed at early 
times. Our simulations further show that the radiative feedback from the first BHs may strongly 
suppress early BH growth, thus constraining models for the formation of supermassive BHs. 
Subject headings: cosmology: observations - galaxies: formation - galaxies: high-redshift - H II regions 
- hydrodynamics - intergalactic medium - black hole: physics 



1. INTRODUCTION 

Radiative feedback from the first generation of stars, 
the so-called Population HI (Pop HI), is a crucial ingre- 
dient in determining how the first protogalaxies assem- 
bled, and in setting the initial conditions for subsequent, 
second-g eneration star formation ([Bark ana & L oeb|2QQl| 
Bromm fc Larson||2QQ4[ [Ciardi &: Terr ara 2QQ5i |Bronim 
et al.||2QQ9p . Ithas been a longstanding goal to un- 
derstand the formation of the first galaxies subject to 
different accompanying feedback effects from the first 
stars, such as ionizing feedback froni individual Pop HI 
stars (iRicotti et al.|2QQ2[ |Abel et al.|2QQ7[ [Yoshida et al 
2007a ), ch emical feedb ack pro duced by a supernova jlSN) 
explo sion (|Wise fc Abel 2008; Greif et al.|2010t|Wise et al.| 
12012), and radiative feedback from accreting black holes 
jBHs) (Kuhlen fc Madau 2005; Alvarez et al. 2009; Park 



RicottT[|201l| . The radiative feedback from accreting^ 



BHs on the formation of the first stars and the assembly 
of the first galaxies is the focus of the current work. 

The first stars are expected to form at a redshift z > 15 
inside dark matter minihalos with masses o f ~ 10*^ Mg, 
Haiman et al.ll996| (Tegmark et al.|l997| [Yoshida et al. 
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Since hydrogen molecules were the only ettec- 
tive low-temperature coolant in the absence of metals, 
the primordial gas remains relatively warm, with typical 
temperatures around 200 K. As a consequence of this in- 
creased thermal pressure, the first stars are expected to 
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be very massive (|Bromm et al.||1999| 2002 |Nakamura fc 


Umemura||2001 


Abel et d.i\[2mi\0mukd.i I 


^ PallallliOOai 


Yoshida et al. 


2006 O'Shea & Norman 2007p. These 



photons, which can photoheat and evacuate the gas re- 
sidin^ within the host minihalos ( Whalen et al. 2004[ 
Alvarez et aL]|2006[ [Johnson et dl|[2^ r ' 



After their short lifetime of ^ 3 Myr, Pop III stars 
end their lives in a SN explosion, or the direct c ollapse 
into a BH. depending on the progenitor mass (Heger 
|et a l.||2003[). If the progenitor's mass lies in the range 
of 15 — 40 Mq, a conventional core-collapse SN occurs, 
or, in the case of rapid rotation, a hypernova. The result- 
ing chemical abundance patterns are in good agreement 
with those obs erved in extremely nietal-poor stars in the 
Galactic halo (iBeers fc Christlieb 2005| Joggerst et al. 
|2010D . 

Xn extreme fate would be a pair-instability super- 
nova (PISN), predicted for progenitor masses between 
140-260 Mo dBarkat et al.|1967||Heger fc Woosley|2002D . 
A consequence of the very large FlISN yields would be 
that even a single PISN could enr ich ^10^ Mr?, of cold, 
dense g^a s u^to lO""^ - 10"^ Z© (iKarlsson et al.||2008 



Greif et al. 2010). Such high metallicities would be m ex- 



cess of any predicted critical value, beyond which a tran- 
sition of star formation mode from massive Pop III stars 
to nor mal Population II (Pop II) would occur ( Omukai 
2000; B romm et al.| |2001a[ [Schneider et al]|2002| : 



For progenitor masses within the range of 40 Mq < 
M < 140 M© or M > 260 M©, the Pop III star wih 
become a massive BH via direct gravitational collapse. 
The ensuing merging of minihalos, provided they con- 
tain cold gas, and the smooth accretion of gas from the 
intergalactic medium (IGM) will feed the central BH. A 
fraction of the accreted mass-energy will then be released 
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Kuhlen fc Madau|2QQ5HWheeler fc jQhnsQrn2Qlli ), 
will ionize the surrounding neutral medium and ne 



as radiative energy, resulting in a miniquasar (iMadau &| 
Rees 2001; Ricotti & Ostriker 2004; Madau et al. 200lf 

whicn 
eat the 
pre-galactic gas. 

Recent simulations have shown that, even in the ab- 
sence of metal-line and dust cooling, the primordial gas 
is able to fragment, leading to a broad initial mass func- 
tion (IM F) including also primordial stars of relatively 
low m ass (Turk et al.'2009l 'Stacy et al."2010"; 'Clark et al." 
2011a b; Greif et al. 2011; Smith et al. 2011; Prieto et al. 



2011[ Stacy et aL|2012p 7which maybe organized in bina- 



ries or stellar systems of still higher multiplicity. Such a 
scenario raises the possibility of the presence of high-mass 
X-ray binaries (HMXB) at z > 10, i.e., binary systems 
composed of a black hole accreting gas from the surface 



of its stellar companion (Mirabel et al. 2011 
HMXBs may provide a much stronger 



Haiman 
eedback 



201 Ij ). 

on the formation of the first galaxies than isolated BHs, 
which have X-ray luminosities that depend crucially on 
the densities and temperatures of the surrounding gas 
jKuhlen fc Madau||2005[ [Alvarez et al.||2009[ [Park fc Ri- 
^otti 201 iTr 

We should note that the assumed presence of HMXBs 
at high redshifts sensitively depends on the Pop III IMF 
and multiplicitly. Recent simulations have provided in- 
triguing hints for the Pop III IMF such that the Pop III 
stars were likely to form in binaries or multiples rather 
than in isolation, but any predictions are subject to 
uncertainties due to the limitations of the sink parti- 
cle method. A sink particle treatment, which has been 
commonly adopted in smoothed particle hydrodynam- 
ics (SPH) simulations to follow the subsequent evolution 
of a protostellar core, allows us to avoid the problem 
of prohibitively s mall numerical timesteps, the so-called 



1995 Bromm et al. 



2002 



^Courant myopia' ([Bate et al. 

Jappsen et al. 200o[ ). 'I'he resulting protostellar llVlI^' is 
attected by the detailed treatment for sink particle merg- 
ers and by the ability of such schemes to self-consistently 
model the intera ction between sink particles and the sur- 
rounding gas ( |Greif et al.||2011 ). It is, therefore, still an 
open question how the early protostellar mass function 
will be mapped into the final IMF of Pop III stars. 

In this paper, we investigate the feedback effects from 
such accreting stellar-mass BHs by focusing on the ques- 
tion: "How does a stellar black hole, a remnant of a 
Pop HI star, infiuence the subsequent star formation and 
in turn the assembly pro cess of the fi rst galaxies?" Our 
work is complementary to Greif et al. ( 2010, ) because our 
cosmological simulations start from tne same initial con- 
ditions. The difference is that we include the feedback 
from accreting BHs while in their work the feedback from 
a PISN was taken into account. We can thus investigate 
how the assembly process of the first galaxies differs in 
response to different assumptions regarding the fate of 
the first stars and the accompanying feedback. In addi- 
tion, we can constrain the resulting gas properties inside 
the centers of such primordial galaxies, thus providing 
us with the initial conditions for second-generation star 
formation. 

It is challe nging to unambiguously define the first 
galaxies (e.g., Bromm fc Yoshida||2011 t). A plausible ap- 
proach is to posit that the first galaxies should be massive 
enough to exhibit a sufficiently deep potential well to re- 



TABLE 1 
Simulation Parameters 



Simulation 


BH feedback 


BH multiplicity 


BHN 


No 




BHS 


Yes 


BH single (isolated) 


BHB 


Yes 


BH binary (HMXB) 



tain gas that is photoheated by stars or BHs. Another 
criterion is that a galaxy should sustain a self-regulated 
star formation mode. Recent cosmological simulations 
suggest that t he fir st galaxies resided in so-called atomic 
cooling halos ( |0h fc Haiman 2002) wit h typical masses 



lO^M^at z 



10 ((W. 



Abelj^2007 



2008 Greffetal 



2008 2010). Whether the presence of an X-ray radiation 



field such as produced by accreting BHs will result in 
positive or negative feedback on the formation of these 
galaxies remains an open key question (|0 h 2001 ; Ricotti] 
eFaI]|2QQl] IGlover fc Brand|[2QQ3l [Machacek et al...2003i 
Ricotti fc O striker 20041 |Kuhlenir Madau|2005[ [Zaroubi 



et al. 2007[ [Ripamonti et aL 2008| Thomas fc Zaroubi 
2008p. A positive feedback would arise if H2 molecule 
formation were promoted due to the elevated fraction of 
free electrons from X-ray photoionization, facilitating the 
cooling and collapse of gas and thus star formation inside 
minihalos. Alternatively, X-ray heating may have sup- 
pressed star formation by increasing the entropy fioor in 
the pre-galactic medium, implying a negative feedback. 

Miniquasar feedback may also be expected to have im- 
portant effects on BH growth. One outstanding question 
is how the supermassive black holes (SMBHs), observed 
by the Sloan Digital Sky Survey (SP SS) at z ~ 6 with 
masses of ~ 10^ Mq (Fan et al.l2006), were able to grow 
within such a short pe riod of time after the Big Bang 
dHaiman fc Loeb||2001[ ). To answer this question, some 
possible scenarios have been suggested, depe nding on the 
initial seed mass dBromm fc Yoshida|[2QTl| . One of the 
main ideas is that the SIVIBH is born as a stellar-mass 
BH, the remnant of a Pop III star with a mass close to 
that of its progenitor, and grows through the accretion 
of surrounding gas and me rgers with B Hs that form in 
neighboring halos ([Madau fc Rees 2001 1 |Li et al.||2007 
Volonteri fc Rees||2006[ [ lanaka fc HaimaS pUpgi r Ho^ 



ever, the photoheatmg and photoionization from a mas- 
sive Pop III star is strong enough to significantly suppress 
further supply of cold gas onto the BH, creating a seri- 
ous early bottleneck to growth |J ohnson fc Bromm 2007; 
Pelupessy et al.|2007| [Alvarez et al., 2009,, Park fc Ricotti] 



20m 

7r"^competing scenario suggests the direct collapse of 
primordial gas into an atomic cooling halo in the absence 
of low-temperature coolants such as H2, or metals and 
dust. Due to the lack of low-temperature coolants, the 
fragmentation of the gas will be avoided, and the ensuing 
isothermal collapse may result in the direct formation of 
a BH without Pop III stellar progenitor. As a result, 
the mass of the seed BH is on the order of 10^ — 10^ M©, 
allowing subsequent growth from a more massive founda- 
tion (^ romm fc Loeb 2003; Koushiappas et al. 2004; Wise 
' " Abel,.2008t ,Regan fc Haehnelt..2009. ,Schleicher etlT 



20101 IShang et al.||2010| [Johnson et al.l|2011| |Latif et al. 
2011"). In this scenario, negative feedback due to photo- 
heating and radiation pressure may still oppose growth 
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Fig. 1. — The emission spectra produced by gas accreting onto a 
100 M© BH with two different accretion rates, M = 10~^ Moyr"^ 
and M = 10'^^ Mqjt-^ . "MCD" (sohd) shows the continuum 
spectra from the multi-color disks and "NT" (dashed) denotes a 
non-thermal component with a power law of /S = 1. The combined 
spectra are normalized to the total luminosity, L = eMc^, where 
e = 0.1 is the radiative efficiency, assuming that each component 
contributes 50% to the total. Note that in our simulations we 
ignore the multicolor disk component and instead assume that the 
radiation emitted by the accreting BH is characterized by a pure 
power law spectrum. 



(Johnson et al. 2011). A second goal of this study is 
to provide an iniproved understanding of whether a stel- 
lar BH can grow to become a SMBH in the presence of 
stellar and BH feedback. 

Here we present the results of cosmological simulations, 
starting from the point where Pop HI stars are first ex- 
pected to form, up to the virialization of a primordial 
galaxy, taking into account both the ionizing radiation 
emitted by Pop HI stars and by miniquasars. We as- 
sume that all Pop HI stars end their lives as massive BHs 
without SN, allowing us to investigate solely the feed- 
back from miniquasars and to compare our results with 
the case where Pop HI stars die as PISNe. The outline 
of this paper is as follows. In Section [2| we describe the 
numerical methodology for the initial setup, the stellar 
radiative feedback, and the feedback from miniquasars. 
In Sectionjsj we show how the primordial galaxy is assem- 
bled with and without BH feedback. Then, in Section [4] 
we discuss the properties of the gas as it is falling into 
the center of the galaxy. We present our conclusions in 
Section [5l 

In the following, all lengths are expressed in physical 
(i.e., not comoving) units unless explicitly stated other- 
wise. 

2. NUMERICAL METHODOLOGY 

In this section, we describe the initial setup of the sim- 
ulations and our methodology for calculating the accre- 
tion rate onto BHs, as well as our treatment of radiative 
feedback from the Pop HI progenitor stars, an isolated 
BH remnant, and from a HMXB. 

2.1. Initial Conditions 



To survey the relevant parameter space, we have car- 
ried out three cosmological simulations. As a refer- 
ence, simulation "BHN" includes stellar radiative feed- 
back from Pop HI stars, whereas the subsequent feedback 
due to BH accretion is not taken into account. The sim- 
ulation "BHS" includes both the feedback from Pop HI 
stars and a single isolated BH. Finally, in simulation 
"BHB" we assume that the Pop HI BH remnant has 
a stellar binary companion, giving rise to a HMXB. We 
provide a summary of the simulations in Table [l] 

For the simulations presented here, we use a parallel 
version of the comb ined hydrodyn amics and tree N-body 
code GADGET2 ( |Springel|[2QQ5t . The code evaluates 
hydro dynamical forces usmg the SPH technique. Our 
simulations are in itialized usin^ a sna pshot from the ear- 
lier simulation of Greif et al. (|2Q1Q[). We start running 
the simulations at z ^ 3U, corresponding to the point 
just before the first Pop HI star in the computational 
box is formed, and terminate them at z ~ 10, when 
the assembly of the first galaxy is expected to be com- 
plete. The original simulation was initialized at z = 100 
in a periodic box of linear size of 1 Mpc (comoving), us- 
ing ACDM cosmological parameters with matter density 
Qrn = 1 — = 0-3, baryon density = 0.04, present- 
day Hubble expansion rate Hq = 70kms~^Mpc~^, spec- 
tral index ns = 1.0, and normalization = 0.9, con- 



sisten t with the WMAP5 measurements ( Komatsu et al. 
2009). Given that we use a higher va lue compared to th e 
current measurement of as = 0.8 (Komatsu et al. 2011|, 
structure formation may be accelerated in our simula- 
tions. However, we compensate to some extent for the 
lack of large-scale power in our simulations due to the 
limited box size, by taking a value larger than the cos- 
mological one. We also stress that in this paper we focus 
on individual halo properties which are not sensitive to 
the variation in the os parameter. 



Using a standard zoom-in technique (see [Greif et al. 
2010), a preliminary run with 64^ particles was hi- 



erarchically refined to generate initial conditions with 
high— mass resolution inside the region destined to col- 
lapse into the first galaxy. Employing three consecutive 
levels of refinement, the mass of DM and gas particles 
in the highest resolution region is moM ^33 Mq and 
^sph ~5 Mq , respectively. The corresponding baryonic 
mass resolution is Mres ^ l-5A/'neigh^sph ^ 400 Mq , 
where Nneigh ^ 50 is the number of neighboring parti- 
cles w ithin the SPH smoothing kernel (Bate & Burkert 
1997). The Jeans mass in primordial gas, where molec- 
ular hydrogen cooling imprints a characteristic density 



of nn = 10 cm and a temperature of 200 K ( [Bromm 
et al. 2002), is thus marginally resolved. 

2.2. Chemistry^ Heating, and Cooling 

We use the same che mistry and cooling network 
(20101), 



Greif et al. 



where all relevant cooling 
mechanisms, such as H and He collisional ionization, 
excitation and recombination cooling, bremsstrahlung, 
inverse Compton cooling, and collisional excitation 
cooling via H2 and HD, are taken into account. 
For H2 cooling, collisions with protons and electrons 
are explicitly included, in addition to the usually 
dominant neutral hydrogen atoms. The code self- 
consistently solves the rate equations for the abundances 
of H,H+,H-,H2,H+,He,He+,He++, and e", as weh as 
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Fig. 2. — Radiative feedback from the first star. We show the 
projected temperature distribution, employing a density-squared 
weighting scheme, of the H II region around a IOOMq Pop III star 
at z = 28, briefiy after its death. The size-*^ of the final H II region, 
~ 5 kpc, greatly exceeds the virial radius of the host halo, ~100 
pc. 

the three deuterium species D, D+, and HD. We can 
thus accommodate the non-equihbrium chemical evolu- 
tion which is ubiquitous in early universe structure for- 
mation. 

2.3. Sink Particle Method 



Emplo ying the sink particle algorithm of Johnson fc 
Bromm|('2007), we convert an SPH particle into a colli- 
sionless sink particle if its hydrogen number density ex- 
ceeds a threshold value of nmax = lO^cm"^. When a 
sink particle forms, gas particles within an accretion ra- 
dius, racc7 are immediately accreted onto the sink. The 
position and velocity of a sink particle are estimated ev- 
ery timestep based on a new mass- weighted position and 
velocity of the accreted gas particles within race = ^res- 
Here, the resolution length of the simulation is defined 
as: 



0.5 



1/3 



Ipc, 



(1) 



where p^ax - ^max^H- 

A sink particle can grow in mass through the further 
accretion of surrounding gas. The criterion for subse- 
quent accretion is that a neighboring SPH particle ap- 
proaches the sink to within race, which we hold constant 
throughout the simulation. Note that our simple pre- 
scription for sink growth is sufficient for the purpose of 
providing a marker for the position of a Pop HI star and 
its BH remnant. When determining the mass growth of 
the BH which is responsible for the X-ray feedback, we 
use a more sophisticated methodology to determine its 
accretion rate (see Section 2.5). 



^ A numerical error in the code for the ionization front tracking 
caused the simulations presented in this work to overestimate the 
radii of the Pop III H II regions by a factor ~ 2, implying final H II 
region radii of ^ 5 kpc instead of ~ 2-3 kpc (compare with, e.g., 
[Greif et al.|2QQ9l their Fi gure 3). However, this overestimate of the 
radii does not affect the conclusions of the present work. This is 
because we are discussing the differences between simulations with 



2.4. Stellar Radiative Feedback 

After a Pop HI star has formed inside its host mini- 
halo, we turn on the radiation field emitted by this star, 
treated as a point source. We propagate the ionization 
front around it to build up a primordi al H II region, usin ^ 
a well-tested ray-tracing algorithm (Greif et al. 2009). 



This scheme solves the ionization front equation m a 
spherical grid by tracking 10^ rays with 500 logarithmi- 
cally spaced radial bins around Pop HI stars. The hydro- 
dynamical effect is taken into account by self-consistently 
coupling the ray-tracing module to the chemical and 
thermal evolution of the gas. 

We assume that a 100 Mq Pop HI star emits blackbody 
radiation with an effective temperature T^ff = 10^-^ K. 
and lumino sity = 10^'^ Lq ( [Bromm et al.||2001b[ 



Schaerer][2Q02 ) . The corresponding production rates tor 



ionizing photons are: A^ion ,Hi/Hei = 9-1 X 10"^^ s ^ and 

^ion,Heii = 4.1 X 10^^ s~^. The growth of the ioniza- 
tion front continues until it reaches its maximum size 
at the end of the Pop HI star's life after t^ = 2.7 Myr. 
However, not all Pop HI stars are endowed with an H II 
region because of the computational expense incurred by 
the ray-tracing. We trace the photons only from those 
stars formed within the Lagrangian volume destined to 
become the first galaxy at z ~ 10, roughly corresponding 
to a 10 kpc radius from the center of this region. 

Finally, we include the transfer of the H2 -dissociating 
Lyman- Werner (LW) photons in the range of 11.2—13.6 
eV, isotropically emitted by a Pop HI star as a 1/r^ field 
without attenuation of the flux with radius, which prop- 
agate far beyond the H II region. The photo-dissociating 
rate is kn^ = 1.1 x IO^Flw? where w is the radiation 



flux integrated over the LW bands (Abel et al. |1997| ). 
We compute Flw in the optically thin limit and do not 
take H2 self-shielding into account (see Greif et al.||2009 



for more details). The neglect of self- shielding is jus- 
tifled due to the low H2 column densities encountered 
in our simulations, of order ^ 10^^ cm~^ inside the BH 



host galaxy ([Draine 
Haiman||2011| r 



Bertoldi 1996 Wolcott-Green & 



2.5. Black Hole Feedback 



In an important difference from I Greif et al. (2010), we 
here assume that all Pop HI stars die via direct collapse 
into BHs, without any preceding SN explosion. For sim- 
plicity, however, we allow only one BH to produce radia- 
tive feedback due to the accretion onto it, taken as the 
flrst BH formed during the simulation. The heat input 
due to this BH feedback into the pre-galactic region is 
active for ~ 350 Myr, while the radiation from individ- 
ual Pop HI stars is turned on for only 2.7 Myr. We next 
describe our treatment for estimating the accretion rate 
onto a BH and the model for the radiation emitted by 
an isolated BH and a HMXB. 

2.5.1. Accretion Rate 

We assume that the BH is embedded in a pressure- 
supported primordial gas cloud, steadily accreting from 

and without black hole feedback, which was correctly implemented, 
and because the radii of the H II regions were overestimated in a 
systematic manner by the same factor in all these simulations. In 
any case, the simulations are consistent with the case including 
feedback from more massive Pop III stars with larger H II regions. 
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Fig. 3. — Gas properties within the H II region depicted in Fig. |2] We show the situation ~1 Myr after the radiation from the Pop III 
star was switched off. From left to right, the panels show gas temperature, hydrogen number density, and electron fraction as a function 
of distance from the central star. The color coding indicates the fraction of mass per pixel in a given phase. The total gas mass, Mtot, 
corresponds to the central 10 kpc. 



it. Tlie correspo nding accretion rate is given by tlie 
Bondi & Hoyle (1944) model, wliere a homogeneous 
medium wliicli is at rest at infinity accretes onto a point 
mass. Tlie Bondi-Hoyle rate can be written as: 



Mbh 



4^(GMBH)Vgas 



sidered tlie case of radiatively efficient accretion. If tfie 
cooling time scale for the liberation of viscously gener- 
ated energy is longer than the accretion time scale, most 
of the energy is advected inward with the accreting gas 
instead of being radiated away ( Nar ayan & Yi|1994|,1995 



^BH 



)3/2 



(2) 



msLeaa oi oemg raaiaLea away i iNa rayan ^ i i|ib>b>4i:||ib>b>0| 
Blandford fc Begelmanll999[ |26u4[ ). A traction ot the en- 
ergy carried by the outflow could ultimately be converted 



where the gas sound speed, Cg, and gas density, Pgas, 
are determined by averaging over the A/'neigh — 50 SPH 
particles closest to the BH. The relative speed of the BH 
with respect to the surrounding gas, ^bh, is negligibly 
small throughout the simulation. Varying A^neigh by an 
order of magnitude does not affect the accretion rates 
when they are estimated with Equation 

By assuming that 10% of the rest-mass energy of the 
accreted matter is released as radiation, we normalize 
the total luminosity according to L = eMBHC^, where 
e = 0.1 is the radiative efficiency and c the speed of 
light. We estimate that the accretion rate Mbh lies be- 
tween 10~^^ M0yr~^ and 10~^ M0yr~^. The former cor- 
responds to the initial situation, where the BH is located 
at the center of an H II region with a temperature of 
^ 10^ K and hydrogen number density of ~ 1 cm~^. 
The upper limit is derived from comparison simulations 
where stellar and BH radiative feedback is neglected. 

It is important to note that the accretion rate used 
in this work is an upper limit, for three reasons. In re- 
ality, when also the feedback from radiation pressure is 
taken into account, the true rates are likely much smaller 
than the nominal Bondi-Hoyle value (Milosavljevic et al. 
2009a|b Park & Ricotti 2012). We estimate, assummg 
a typical accretion rate of 10~^ M0yr~^, that the accel- 
eration due to radiation pressure from Thomson scat- 
tering and photoionization in the vicinity of the BH, at 
r = 1 pc, is 30 times larger than the gravi tational accel- 



eratio n, ttgrav = GM^n/r^ (e.g., Eq. 10 in [Johnson et al. 
2011[ ) . This indicates that radiation pressure, mamly due 
to photoionization, is important near the BH, acting to 
limit growth. The other reason is that we have only con- 



y carriea by the outflow could ultimately 
to radiation, but the resulting spectrum would certainly 
be very different than that of the thin disk assumed here. 
A third reason is that the vorticity of the turbulent gas 
has not been proper ly resolved here. Studyi ng its effect 
on Bondi accretion, Krumholz et al. (2005) found that 
even a small amount ot vorticity is able to significantly 
reduce the accretion onto the BH. This suppression could 
be important in an atomic cooling halo, where the radial, 
cold flow along the filaments is converted into turbulent 
motion, thus generating vorticity in the central region of 
the galaxy. 

In addition to the Bondi-Hoyle prescription, there is 
another way to estimate BH growth, based on Mgink 
calculated with the method in Section |2.3[ where the 
BH mass grows by accreting all gas particles within 
race- The resulting rate is much higher than the Bondi- 
Hoyle value, roughly by a factor of 10. We adopt 
Mbh rather than Mgink as an estimator for BH growth, 
given that the sink accretion radius is about one or- 
der of magnitude larger than the Bondi-Hoyle radius, 
tb = {/j^ttiylG Mbh) / ^ where /i is the mean molecu- 
lar weight and /cb is the Boltzmann constant, or, tb ~ 
O.Spc (Mbh/IO^Mq) (T/102K)"\ Moreover, only a 
fraction of the infalling gas will likely reach the vicinity 
of the BH. Furthermore, we can thus track the change in 
BH mass in a smooth manner, and do not have to con- 
tend with discreteness effects owing to the limited mass 
resolution afforded by SPH. 

2.5.2. Non- Thermal Radiation 

In the local universe, the emergent spectra from ac- 
creting BHs are typically modelled as a combination 
of a power-law component, F^, ex z/^, describing non- 
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thermal (NT) synchrotron radiation, and a multi-color 
disk (MCD), resulti ng in a soft^hermal continuum (e.g., 
Mitsuda et al.|1984 ^. In Figure [l] we show typical spectra 
to illustrate this model where a soft thermal component 
is fitted by a MCD blackbody with kT^^ - 150 eV (Miller] 
|et al. 2003), normalized assuming two representative ac- 
cretion rates. Here, Tin is the temperature of the inner 
disk region which is related to the BH mass according 

It is currently not known whether ac- 



to 



-1/4 
BH 



creting BHs at high redshift would behave in a similar 
manner. 

Conservatively, we here assume that there is no signif- 
icant difference in the BH emission physics over cosmic 
time, and use the same spectra. However, we ignore the 
MCD component of the emitted radiation in our simula- 
tions. Such a thermal contribution would determine the 
radiative feedback in the immediate "near zone" from the 
central BH, but is unimportant on larger scales, where 
the NT component dominates. We estimate the size 
of the near zone with a standard Stromgren analysis. 
Assuming that the ionizing photons originate from the 
MCD component, and employing a typical accretion lu- 
minosity of Mbh = 10~^ Mq yr~^ and a typical density 
of nH = 10 cm~^, we compute an equilibrium ionization 
front (I-front) radius of < 5pc, which is comparable to 
the spatial resolution of the simulat ions. 

We follow [Kuhlen fc Madau] ( |2005| to model the prop- 
agation of high-energy photons where the isotropic radi- 
ation field, (X 1/r^, is a function of distance from the BH 
only. The coefficient for the ionization rate, /cion, can be 
written as 



r 



(3) 
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Fig. 4. — Distances between newly formed Pop III stars that are 
accompanied by H II regions and the BH as a function of redshift 
from our three simulations. Note that the BHN, BHS and BHB 
simulations end at different redshifts, corresponding to z ~ 11.6, 
z ~ 10.4, and z ~ 15.8, respectively (indicated by the vertical dot- 
ted lines). A total of ~ 50 Pop HI stars form within 10 kpc around 
the central BH in the range of redshifts z = 20 — 10 in simulations 
BHN and BHS. In the BHB run, the positive feedback is clearly 
evident far away from the source. This is due to the gas collapse 
into distant minihalos, facilitated via H2 cooling promoted by the 
strong X-ray emission from the HMXB; locally, on the other hand, 
star formation is suppressed due to the strong negative feedback 
from the binary source. We also show the virial radius of the DM 
halo hosting the active BH or the HMXB((ia5/ie(i line). The X- 
ray photoionization heating from a stellar-mass BH is also able to 
quench further star formation in the host halo at all times before 
the halo enters the atomic cooling phase. The top axis shows the 
time elapsed since the formation of the black hole at time to- 



where z^nt = 0.2keV//i is the frequency where the 
non-thermal component starts to dominate and z^max = 
10keV//i is the upper cutoff frequency. Our choice of 
energy range , 0.2 — lOkeV, is motivated by the study of 



Miller et al. 



j2QQ3[ ), where the spectra of select ULXs 
associated with lUOM© black holes were fitted with a 
power law in this same range. For the hydrogen and he- 
lium photo-ionization cross sec tions, a^, we use the stan- 
dard expressions (e.g., Barkan a fc Loeb|2QQl Osterbrock 
[fc Fer land 2006). 

reconsider the secondary ionization effects from ener- 
getic electrons released through the absorption of an X- 
ray ph oton, we adopt the fitting formulae given by|Shull 



van Steenberg (1985) (see also Valdes fc Ferrara||2008 
Furlahetto & Stoever 2010). They computed the Irac- 



tions of the initial electron energy going into secondary 
ionizations of HI, secondary ionizations of Hel, and into 
heating the surrounding gas. Secondary ionizat ions of 
Hell are not important ( |Shull fc van Steenberg| [T985), 
and are thus not included in the overall energy budget. 
The effective ionization rates can thus be written 



where 



fceff[HI] =kio„[HI]+k,ec[HI], 



fc,ec[HI] =fH (rHI + — FHel^ ^ 



nm / 13.6 eV 



(4) 



(5) 



and 



where 



/ceff [Hel] = kion[HeI] + k^ecpel]. 



/Csec[HeI] = fHe ( ThcI 



hhi j 

riHel 



HI 



1 



24.6 eV 



(6) 



(7) 



Here F, the rate at which the excess energy is released 
due to the first photoionization of the gas by X-rays, is 
given by 



, / F^cr^ 1 — dv, 



(8) 



where nn is the number density of the respective union- 
ized species, i^min the ionization threshold frequency, 
^^min=13.6 eV, /ii^min=24.6 cV, and /iz/niin=54.4 eV for 
hydrogen, neutral helium, and singly ionized helium, re- 
spectively. 

The fractions going into secondary ionizations, /h and 
/He, are a function of the hydrogen ionization fraction, 
^ion = ^H+/^H, and can be expressed as 



/h= 0.3908(1 -xl^.r')'-'''' (9) 

/He = 0.0554(1 (10) 
and the photo- heating rates are then H = r(l — /n/He)- 
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Fig. 5. — Top panel: Comoving stellar mass density vs. redshift. 
We compare the results from our three simulations, the BHN com- 
parison case {solid line), the BHS feedback run (dashed line), and 
the BHB run (dot-dashed line). It is evident that the positive 
feedback from the BH binary, acting on cosmological scales, has 
boosted star formation at the highest redshifts. Bottom panel: Co- 
moving star formation rate density vs. redshift. The lines refer to 
the same simulations as above. For comparison, we overplot the an- 
alytic fitting formula, derived by Hernquist &; Springel (2003), for 
higher-mass systems where atomic cooling is active [dotted line), 
and that of Yoshida etaL. C2003) for Pop III star formation in 
minihalos via II2 coolmg [dashed triple-dotted line). 



The flux from the non-thermal component is given by 



^0 



-/3 



/3 = 1, 



(11) 



in the range of 0.2—10 keV, implying a mean hydrogen- 
ionizing photon energy of 390 eV. Here F^^q, the flux at 
^0 = ^NT, is determined at every timestep, based on the 
normalization: 



cMbhc^ 



(12) 



where, L^, = A.nr'^Fy is the specific luminosity. We note 
that since extinction along the line of sight is not taken 
into account, the estimated ionization and heating rates 
are likely upper limits. Recall that the accretion rate, 
Mbh, is updated at every timestep, thus allowing our al- 
gorithm to reflect the changing conditions in the vicinity 
of the BH. 

The resulting photo-ionization and photo-heating rates 
are a function of distance from the BH, the instantaneous 
accretion rate, a function of the hydrogen ionization frac- 
tion, and can be conveniently written as 



pc 



Ml 



BH 



10-6 Mo /yr 



(13) 



K = [1.96, 2.48, 0.49] x 10"^^ s"\ (14) 



and 



/Hej , 



where 



(15) 



H = [7.81, 9.43, 1.63] x 10"^^ ergs"^ (16) 



for neutral hydrogen, neutral helium, and singly-ionized 
helium, respectively. The corresponding number densi- 
ties of each species j are given by nj . 

2.5.3. HMXB Emission 

Recent improved simulations of Pop HI star formation 
suggest that the primordial gas could fragment into two 
or more distinct cores, possibly resultin g in a massive bi- 
nary, or higher-m ultiple ste llar syste m fTurk et al. 2009 
Stacy et al.||2QTQl [20 12; Cl ark et al.l, 2011a b; |(jrell et af 



20TTr Smith et al, 



ture provides anoi 



et al. 


2011 


Haiman 2011) 


bounc 


1 to the slightly less 



pOlli [Prieto et al.||2011| ). This pic- 
ner possible X-ray source connected 
to the first stars, where the primary in a binary system 
leaves a BH behind at the end of its short life (Mirabel| 

The remnant could remain 
assive secondary, where nu- 
clear burning is still going on. The material lost by the 
secondary in a wind or by Roche lobe overflow will be 
dumped onto the companion BH, releasing copious X- 
ray emission. The duration of such a HMXB is limited 
by the main-sequence lifetime of the donor, and is of the 
order of 10^ yr. 

It has been suggested that there is a correlation be- 
tween the number of HMXBs and ultra-luminous X-ray 
sources (ULXs), which can be explained as a subpop- 
ulation of the former, and lo w-metallicity environment 
(iMaiid et a l. 2004; Dray 200^ ISoria||2007| [Mapelh et al. 
z009[ Linden et al.| |2010|) . We can thus expect that 
the number densities ol HMXBs and ULXs at high red- 
shifts, when the universe was chemically pristi ne, were 
higher than wh at is observed locally ( Mirabel et al.|2011 



Haiman 



2011). Such local surveys show that HMXBs 



and ULXs have luminosities of Lx ~ 10^^ erg s ^ and 



Lx > 10^^ 



erg s 



respectively, which are clos e to the 



Eddington value, or mildly super-Eddington (Grimm 
et al. 2003). Consequently, we assume that the lumi- 
nosity of a HMXB at high redshifts is close to the Ed- 
dington limit, LEdd = 47rGMBH^Hc/crT, where mn, ctt, 
and c denote the proton mass, Thomson cross section, 
and the speed of light, respectively. We further assume 
that the spectrum exhibits the same power-law behavior 
as in Equation ( 11 ). The X-ray luminosity in the HMXB 
scenario is then oetermined from 



LEdd = 1.38 X 10^^ erg s"^ 



Ml 



BH 



100 Mo 

and we assume Mbh = 100 M©. 

3. FIRST GALAXY ASSEMBLY 
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In this section we present the results obtained from 
our simulations. In Section |3.1[ we briefly discuss the 
properties of the H II regions around massive Pop III 
stars in order to check for agreement between our work 
and previous studies. We proceed in Section 3.2 to de- 
rive the star formation rate in the Lagrangian volume 
of the emerging primordial galaxy. We then discuss in 
Section [33| how this galaxy is assembled, contrasting the 
simulations of feedback from a single isolated accreting 
black hole (BHS) and from a HMXB (BHB) described in 
Sec.[2]with a simulation without BH feedback (BHN). We 
conclude with a description of the accompanying black 
hole growth. 

Note that the BHN, BHS and BHB simulations end at 
different redshifts, corresponding to z ^ 11.6^ z ^ 10.4, 
and z ^ 15.8, respectively. The reason that we termi- 
nated the BHN run earlier than the BHS case is mainly 
due to computational expense, but the galaxy at z 11.6 
already meets the criterion for first galaxies and the re- 
sults give us enough information to compare with those 
of the BHS run. For the BHB case, we have performed 
the simulation for much longer than any realistic lifetime 
of HMXBs. For example, the main-sequence lifetime of a 



donor star of 30 M^ , close to the predicted ty 



mass ( Clark et al. 2011b: Hosokawa et al. 
^^Im\] wouM be ^ 20Myr. 



3ical Pop HI 



2011 Stacy 



et al. 2012p , would be ^ 201VLyr. 'I'hus, our work can 
be considered to represent the case where at least one 
HMXB exists continuously over 150 Myr, somewhere in 
the Lagrangian volume of the emerging protogalaxy. 

3.1. H II Regions around the First Stars 

The first Pop III star to appear in our cosmological box 
forms at a redshift z = 28 inside a 5 x 10^ Mq minihalo. 
Since the Pop HI star is assumed to be very massive, 
with 100 M0, it produces copious amount of ioniz- 

ing photons which carve out an extended H II region into 
the surrounding IGM. The gas within the H II region is 
photo-heated to temperatures, T ~ 3 x 10^ K, in excess 
of the virial temperature of the minihalo, thus triggering 
a hydrodynamic outflow. In Figure [2j we show the pro- 
jected temperature approximately 1 Myr after the radia- 
tion from the central star was switched off. The terminal 
size of the H II region is about ~5 kpc^, such that the 
I-front extends much farther out than the virial radius 
of the minihalo, 

' vir 

100 pc. Note that the anisotropic 
morphology of the H II region reflects the inhomogeneous 
density distribution in the neighboring IGM. We describe 
the physical properties of the photo-ionized gas in Fig- 
ure [3j again ~ 1 Myr after the Pop HI star died. At this 
time, the gas towards the center has already begun to 
cool and recombine, leading to a relic H II region. 

The expansion of the I-front begins with a very short 
lived R-type (supersonic) phase, followed by a D-type 
(subsonic) one, where the I-front is trapped behind a 
hydrodynamical shock wave (e.g 



Alvarez et al. 



2006 



yr, the I- front is ab 



Johnson et al. 



Whalen et al. 2004 



2007). After rZlQ^ 



e to break out, now again supersoni- 



cally racing ahead of the shock into the low-density IGM. 
The resulting density profile, together with the location 
of the shock at ~200pc, and a corresponding velocity 
of '^sh ~ 25 — 35 km s~^, are in good agreement with 
the analytical champ ague- flow solution ( |Shu et al.|2002, ). 
We note that contrary to the minihalos of ^ 10*^ M©, 
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Fig. 6. — Assembly of the first galaxy. Shown is the redshift 
evolution of DM virial mass (black solid line) of the halo which 
will host the first galaxy at z ~ 10, as well as of the total gas mass 
(red) and the cold gas mass (cyan) within the halo for simulations 
BHN (solid lines), BHS (dotted lines), and BHB (dot-dashed line). 
The black dashed line represents the critical mass required for the 
onset of atomic cooling in the halo at a given redshift. There is 
no cold gas within the halo in simulation BHB owing to the strong 
heating from the HMXB. 

common at high redshifts, the extent of the shock ra- 
dius is likely to decrease in rare higher-mass halos, with 
> 10^ M0, due to their deeper potential wells and the ef- 
fect of infalling gas, possibly trapping the I-front within 
the host halo, resulting in a compact H II region (Ki- 
tayama et aT||2004[ [Yoshida et "ar||2007aD . 

At the end of the life of the Fop 111 star after ^2.7 
Myr, when it leaves a massive BH remnant behind, the 
surrounding gas has been evacuated from the minihalo, 
thus preventing the relic BH from accreting any cold gas 
for an extended period of time. After the radiation from 
the central star shuts down, the relic H II region starts 
cooling partly due to adiabatic expansion cooling, and 
partly due to atomic hydrogen line cooling, facilitated 
by collisional excitations from the enhanced abundance 
of free electrons. The elevated electron fraction catalyzes 
the formation of H2 and HD molecules in the relic H II 
region, thus providing additiona l molecular cooling (| Ri-| 
cotti et al."2001j [Johnson et al.||2007 ). Roughly 10 
after the radiation from the Fop 111 star turns off, the 
temperature decreases from 10^ K to 10^ K and the H2 
abundance increases from 10~^ to 10~^. While H2 and 
HD molecules are susceptible to photo-dissociation from 
Lyman- Werner (LW) radiation produced in neighboring 
star- forming halos, the relatively short time required for 
the re- formation of H2 and HD, compared to the average 
time between the formation of Pop HI stars, allows their 
abundances to remain high. 

3.2. Star Formation 

Figure [4] shows the distance between newly formed 
Pop HI stars and the central BH in the range of red- 
shifts z = 10 — 20 from three simulations. A total of 
~ 50 Pop HI stars have been formed in the BHN and 
BHS cases, accompanied by individual H II regions, ac- 
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Fig. 7. — Gas properties in the vicinity of the accreting BH in simulation BHS. Shown are the density-weighted temperature (top row) 
and the projected hydrogen number density (bottom row) of the gas within cubical slices of linear size 100 kpc (comoving). The small 
insets provide zooms into the central 10 kpc (comoving). The situation is depicted at three different times, At= 0, 100, and 300 Myr after 
the BH has been formed at the center of the minihalo, corresponding to z = 27.9, 17.7, and 11.0, respectively. The location of the BH is 
marked by a black circle at the center of each plot. 



cording to our criterion that stellar radiative feedback is 
taken into account only for those stars formed within a 10 
kpc radius from the center of the emerging protogalaxy. 
We find that for over 250 Myr after the seed BH formed 
there has been no further star formation within the host 
halo in simulation BHS. This is because it takes time for 
the gas expelled by the BH progenitor star to be reincor- 
porated into the halo, and the modest feedback from the 
BH prevents the gas from cooling. In simulation BHN, 
on the other hand, stars inside the host halo continue 
to form already much earlier, demonstrating that X-ray 
heating from the accreting BH in simulation BHS implies 
a strong local negative feedback. 

Note the intense burst of star formation taking place 
within the virial radius of the protogalactic halo at z = 
11.5 — 10.5 in the BHS. This starburst is fueled by a 
massive infall of gas. This gas is able to cool not only via 
molecular hydrogen cooling but also via atomic hydrogen 
line cooling. A fraction of the gas mass moving toward 
the central BH is thus consumed by star formation rather 
than being accreted onto the BH. 

In Figure js] we show the evolution of the stellar mass 
density and the star formation rate density (SFRD). The 
comoving stellar mass density is calculated using the stel- 
lar masses that have been formed within the 10 kpc ra- 
dius from the center of the emerging protogalaxy at a 
given redshift. For the SFRD, we take the time deriva- 
tive of the stellar mass density at the redshift when 
the Pop HI stars formed. For comparis on, we over- 



Springel 2003), for higher-mass halos where atomic hy- 
drogen cooling is dominant, and for minihalos, assuming 
that on e 100 Mc?) Pop III st ar forms per system via H2 
cooling ( Yoshida et al.| 2003). 

The estimated star formation rate densities from our 
simulations lie between the two analytical fits. Evidently, 
at z = 20 the SFRD in the BHB run is higher than 
in the BHN and BHS cases, by a factor of about 3. 
This is a consequence of the positive feedback, where 
gas collapse into distant minihalos is facilitated via H2 
cooling promoted by the strong X-ray emission from the 
HMXB. As density fluctuations grow and star formation 
is enhanced, the local negative radiative feedback from 
Pop III stars begins to dominate, mitigating the positive 
HMXB feedback effect. We have not considered cooling 
due to metals and dust, both produced in SN explosions. 
At low densities, n < 500 cm~^, and at low metallici- 
ties, Z < 10~^ZfT^, however, H2 coolin g is expected to 



dominate over metal-line cooling ( jJappsen et al. 2007 
,2009a,bJ . 

3.3. Mass Growth 

The characteristic property to distinguish the first 
galaxies from lower mass minihalos is their ability to cool 
through atomic hydrogen line emission, depending on the 
virial temperature of the hosting halo with a virial mass 
Mvir: 



plot SFRDs using analytic fitting formulae ( Hernquist 
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Fig. 8. — BH growth with and without feedback. Shown is the 
redshift evolution of the BH accretion rate, the temperature and 
density of the gas in the immediate vicinity of the BH, as weh as the 
resulting BH masses for simulations BHN (black) and BHS (red). 
In the top panel, we also indicate the corresponding Eddington- 
limited accretion rates for the two cases (dashed lines). 

Above Tvir ^ 10^ K, the gas within the halos is able to 
cool mainly via atomic hydrogen. In Figure [6j we show 
the evolution of the virial mass for the most massive halo, 
which will host the first galaxy at z ~ 10, as well as of 
the total and cold gas masses for the three simulations 
presented in this work. The virial mass of the DM halo 
is estimated as the mass within a sphere with average 
DM density p ~ 200 po{z)^ where po is the mean cosmic 
density at a given redshift. Cold gas is identified by the 
condition that the gas temperature is less than half of 
the halo virial temperature, T < 0.5 Tvir(^). We also 
indicate the critical mass required for the onset of atomic 
hydrogen cooling at a given redshift. 

We find that at z > 18, the halo is dominated by hot 
gas, exceeding the amount of cold gas by an order of 
magnitude. As time passes on, the cold gas mass in- 
creases, eventually accounting for > 80% of the total gas 
mass in simulations BHN and BHS. This trend can be 



understood by the vulnerability of the halo gas to stellar 
radiative feedback. The corresponding evacuation of gas 
from the halo is very strong at high redshifts, z > 18, be- 
cause the halo potential wells were not yet deep enough 
to retain photo- heated gas. We note that the sharp dips 
in the amount of cold gas are due to star formation in- 
side the halo itself (at z < 12), as well as in neighboring 
halos that are sufficiently close, within < 5 kpc. The ion- 
izing feedback from the accreting BH starts operating at 
z < 18, above which the effect is too weak as a result of 
the initially very low BH accretion rates. 

While the total amount of gas is not sensitive to the 
BH feedback, the reduction in cold gas mass by a factor 
of ~ 5 indicates that the additional heating from this 
feedback, on top of the stellar feedback, has a significant 
impact on the gas in the center of the forming galaxy. As 
the halo grows further via smooth accretion and mergers 
with minihalos, however, at z ~ 13, both the total gas 
mass and mass of cold gas are no longer sensitive to the 
BH radiative feedback. At z ^ 11.5, the condition for an 
atomic cooling halo is satisfied. 

For simulation BHB, the heating from the HMXB 
is so strong that all gas particles have temperatures 
T > 0.5 Tvir (2^), over the entire range of simulated red- 
shifts z > 15, and the total gas mass is reduced by an 
order of magnitude by photo-evaporation, as is clearly ev- 
ident in the dot-dashed lines Figure [6) This result would 
imply that if an HMXB existed within a minihalo at high 
redshifts, it would take significantly longer for the halo 
to reassemble the lost gas, and to eventually evolve into 
a primordial galaxy. 

In Figure [7| we show the evolving conditions of the gas 
in the vicinity of the accreting BH in simulation BHS, 
while the first galaxy is assembled. Initially, the BH 
is located at the center of the H II region produced by 
the progenitor star. As filamentary structure develops, 
gas infall and mergers with minihalos frequently occur 
along the filaments (see the middle panel of Fig.ul. Cold 
accretion along the filaments, which are dense enough to 
allow molecule re-formation, in turn leading to enhanced 
cooling, efficiently delivers cold gas into the center of the 
halo. Roughly 300 Myr after the seed BH formed, the 
conditions for the onset of atomic hydrogen cooling are 
met. The virial radius of the halo at this time is Tvir '^l 
kpc. 

3.4. Black Hole Growth 

Figure |8] shows the evolution of the accretion rate onto 
the BH, the temperature and density of the neighboring 
gas, as well as the BH mass for simulations BHN and 
BHS. For the temperature and density, we average over 
the 50 closest SPH particles. The BH remnant, originat- 
ing in the death of a Pop HI star, initially resides within 
a medium of high temperature, ~ 10^ K, and very low 
density, nn ~ 0.1 cm~^, preventing the gas from cooling 
and hence accreting onto the BH. At this stage, the accre- 
tion rate is only Mbh ~ 10~^^ M0yr~^, and the corre- 
sponding accretion luminosity, assuming a 10% radiative 
efficiency, is Lace = cMbhc^ ~ 5.7 x 10^^ ergs~^. This 
value is seven orders of magnitude below the Eddington 
luminosity, and BH growth and feedback are negligible 
in the beginning. As the potential well of the host halo 
becomes deeper, however, the amount of infalling gas in- 
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Fig. 9. — Gas properties within the virial radius of the first galaxy halo. We compare the situation in simulations BHN a^t z = 11.6 
(left-hand) and BHS at z = 10.4 (right-hand) , corresponding, respectively, to the end of each of the two simulations. From top to bottom, 
each panel shows the gas temperature, electron fraction, H2 abundance, and HD abundance, respectively. The feedback from the isolated 
BH significantly modifies the thermodynamic properties of the central gas. The color coding in these diagrams indicates the fraction of 
mass in a given phase. 



creases with time, boosting the accretion rate, especiahy 
at z ~ 17 — 19, as can be seen in Figure|8] Consequently, 
the accretion luminosity also becomes large enough to in- 
fluence the surrounding gas, keeping it at a temperature 
that is an order of magnitude higher than in the BHN 
comparison simulation. 

The accretion rate in the BHN simulation is already 
comparable to the Eddington value, depicted as dashed 
lines in Figure [s] (top panel), at z ^ 15.5, while it is still 
an order of magnitude lower in the BHS case. Occasion- 
ally, star formation takes place very close to the BH, e.g., 
~ 1 kpc away at 2; ~ 14 (see Fig. [4|. The radiative feed- 
back from this event acts to compound the heating effect 
from the BH accretion, thus rendering the removal of 
gas out of the shallow potential well more effective. Such 
gas evacuation at z ~ 14 due to both feedback effects is 
clearly seen in Figure [6j 

The combined stellar and BH radiative feedback results 
in an accretion rate that is on average 4 orders of magni- 
tude below the Eddington value at z = 14 — 13. Even 300 
Myr after the formation of the BH, its mass has increased 
by only 1.5 % in the BHS simulation, whereas there is a 
two orders of magnitude growth in the BHN case. This 
indicates that the feedback from a stellar-mass BH is suf- 
ficiently strong to prevent significant growth, suggesting 
a very important constraint on SMBH formation scenar- 
ios. We infer that the radiative feedback from an accret- 
ing BH might be partly responsible for the low density 



of quasars at redshifts z 
growth. 



6, by suppressing early BH 



4. PROTOGALACTIC GAS PROPERTIES 

In this section, we discuss the properties of the gas as 
it falls into the center of the emerging galaxy at z ~ 10. 
We compare the situation for our different assumptions 
on the feedback acting during the preceding assembly, 
including a discussion of the complementary PISN case. 
This allows us to constrain the initial conditions for the 
subsequent starburst inside the first galaxies, a crucial 
input for predicting their observational signature. 

4.1. Isolated Black Hole Case 

In Figure [9) we compare the temperature, electron frac- 
tion, H2 fraction, and HD fraction of the primordial gas 
within the virial radius of the first galaxy halo as a func- 
tion of number density for simulations BHN and BHS. 
We see that at low number densities, nn = 10~^ — 1 
cm~^, corresponding to the outskirts of the host halo, 
for both the BHN and BHS simulation the gas tempera- 
ture is high, T ~ 10^ — 10^ K, and the electron fraction 
is enhanced, f^- = 10~^, due to the accumulated relic 
H II regions produced by previous star formation inside 
and outside of the halo, and due to heating in accretion 
shocks. 

In the BHS run, additional heating from the BH al- 
lows the gas within the halo to maintain a relatively high 
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10^ 10^ 10^ 10^ 10^ 

t [yr] 

Fig. 10. — Chemical evolution of the central gas cloud. We calcu- 
late an idealized one-zone model, tracing the history of a represen- 
tative fluid element in the center of the emerging galaxy, starting 
from T > 10^ K and subsequently cooling down. We employ the 
same chemistry module as used in the three-dimensional simula- 
tions. Top panel: evolution of the particle abundances. Bottom 
panel: Characteristic timescales. In the absence of feedback from 
a BH (simulation BHN; solid curves), the H2 abundance converges 
to the asymptotic value of ~ 10~^, reached when both the 
timescale for the evolution of electron abundance and the cool- 
ing time become shorter than the H2 formation and dissociation 
timescales (Oh &; Haiman 2002). With BH feedback present (BHS; 
dashed curves), the electron fraction is enhanced due to the X-ray 
photoionization, and the H2 fraction is enhanced, reaching up to 
> 10~^. We note that the additional photo-heating and photo- 
ionization due to the radiation from a BH increase the cooling 
time and the timescale for the evolution of the electron fraction, 
leading to the peaks clearly seen in the bottom panel. 

temperature over the entire range of densities. There is a 
substantial difference seen in the electron fraction: while 
varying over the large range of /g- = 10~^ — 1 in the BHN 
simulation, in the BHS case most gas within the virial 
radius has electron fractions above > 10~^, and there is 
no gas present with fractions < 10~^. In response to the 
substantially increased electron fraction, the formation of 
H2 and HD is catalyzed, leading to molecule abundances 
that are higher by two orders of magnitude compared to 
the BHN run. 

In the absence of feedback from a BH in simulation 
BHN, the molecular hydrogen fraction converges to the 
well-known asymptotic value of . fa^ ^ 1Q~^. Thi s be- 
havior reflects a freeze-out process ( IOhfcHaimaji|2 QQ2), 
where the molecular hydrogen abundance no longer 
evolves once the recombination and cooling timescales 
become less than the H2 formation and dissociation 
timescales, trec^^cooi <^ ^form, ^diss- Here, the recombi- 
nation time is tree ^ l/(<^B^/e-)5 and the H2 formation 




-6.5 -5.9 -5.2 -4.6 -4.0 



Fig. 11. — Gas properties in the center of the emerging first 
galaxy. We show temperature, chemical abundances, and number 
density vs. radial distance from the galactic center for the same 
simulation BHS at z = 10.4 as in the right-hand panel of Figure lo] 
The radiative feedback from the accreting BH is responsible for the 
presence of two distinct branches within the central 100 pc. 

and dissociation timescales are tform ^ ^H2/^H2,form and 
^diss ~ ^H2/%2,diss, respectively. 

To illustrate the basic difference between the BHN and 
BHS simulations, we calculate a simple one-zone model, 
using the same chemistry module as in the simulations, 
but focusing on a single, representative fluid element as it 
falls into the center of the emerging galaxy. In the model, 
the radiation from the BH is imposed by a 1/r^ radiation 
field, adopting a typical accretion rate of lO~^M0yr~^, 
and we assume that the fluid element has a fixed density 
of n = 10 cm~^, placed r = 100 pc away from the BH 
at z = 11.6. Our fluid element is initially hot, T = 10^ 
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Fig. 12. — History of a representative SPH particle in the vicinity 
of the central BH. We follow the evolution for 60 Myr, starting at 
a distance of 100 pc away from the BH, all the way to the center 
of the emerging first galaxy halo. The spikes in the temperature 
and electron fraction are due to individual star formation events 
occurring inside or outside of the growing halo. As the gas par- 
ticle moves to the center of the halo, feedback from the accreting 
BH begins to operate, promoting H2 and HD molecule formation. 
Very close to the BH, however, the heating is sufficiently strong to 
destroy the molecules, generating the two distinct branches for r < 
50 pc seen in Fig. [iT] 

K, and highly ionized, typical for a location inside a relic 
H II region. The results of the one-zone model are shown 
in the top panel of Figure [lO] To guide the interpreta- 
tion of these results, we also show, in the bottom panel 
of Figure [lOj the cooling time and the timescales for for- 
mation and dissociation of molecular hydrogen, and for 
the evolution of the electron fraction, te = Ue/he. 

In the case without BH feedback (solid lines), the H2 
formation and dissociation timescales are significantly 
longer than both the timescales for the evolution of the 
electron fraction and the cooling rate once the temper- 
ature has dropped below a few 10^ K, establishing the 
freeze-out value of ~ 10~^. In the presence of BH 
feedback (dashed lines), on the other hand, the addi- 
tional ionization from the BH X-ray emission signifi- 
cantly increases the fraction of free electrons at low tem- 
peratures. The enhanced electron fraction promotes the 
formation of H~, and this, at low temperatures, in turn 



implies an increase in the H2 fraction. We note that sim- 
ilarly high H2 fractions have been found in the related 
but different context of the non-equilibrium chemistry 
and cooling of primordial gas behind structure forma- 
tion s hocks under external irradiation ( [Shapiro fc Kang 
19871) 



The overall distribution of the HD abundance is sim- 
ilar to that of H2. In the BHS run, the HD frac- 
tion increases to /hd ^ 10~^ which exceeds the critical 
level needed for efficient cooling to the cosmic microwave 
background (CMB) temperature in local thermodynamic 
equilibrium, .faD.crit ^ 1Q~ ^, by two orders of magnitude 



( Johnson fc Bronini| 2006 ) , as can be seen in Figure!^ 
iNevertheless, HD cooling does not succeed in tying the 
temperature to the CMB because of the continuous heat- 
ing from the BH. 

We plot the properties of the gas in simulation BHS in 
Figure [11] but this time as a function of distance from 
the central BH. The most noticeable feature here is the 
presence of the two distinct branches towards the center 
r < 100 pc which is explained by the BH feedback. To 
better understand the BH feedback effect, in Figure 12 
we follow the evolution of one representative SPH par- 
ticle for ~ 60 Myr. Initially, the tracked particle is lo- 
cated at 100 pc away from the central BH, where BH 
feedback is negligible. The sudden jumps in both tem- 
perature and electron fraction indicate a number of star 
formation events inside and outside of the halo. Soon 
after such events, the gas begins to cool back to a few 
~ 10^ K, as is typical for relic H II regions. 

As the gas particle moves within 100 pc from the BH, 
however, the continuous heating from the accreting BH 
starts to operate, resulting in temperatures of ~ 10^ K 
and an electron fraction of /g- ~ 0.1. As a consequence, 
H2 molecule formation is catalyzed, leading to fractions 
of /hs ^ 10~^. Finally, when the gas particle is located 
too close to the BH, r < 10 pc, the strong X-ray ffux heats 
up the gas again to ~ 10^ K, leading to the dissociation of 
molecules and producing the two distinct branches within 
the central 100 pc, clearly seen in Figure [TT] We point 
out that the dense clumps around 700 — 800 pc from the 
BH in Figure pT] correspond to the middle branch in the 
H2 and HD molecule fractions, evident in the right panel 
of Figure [9| 

4.2. Black Hole Binary Case 

In simulation BHB, where the continuous feedback 
from an HMXB is included, a dramatic change is seen in 
the gas properties surrounding the central HMXB. Fig- 
ure [13] shows the comparison between simulations BHN 
and BHB at z = 15, plotted in the same fashion as in Fig- 
ure [11] but now extending to 30 kpc from the HMXB to 
illustrate the effect of feedback not only on the gas within 
the host halo, but also on the distant IGM. We carefully 
select the snapshots to guarantee that we are seeing the 
feedback effect from the HMXB, without confusion from 
Pop HI stellar radiation. In the BHB run, the gas tem- 
perature within the central 1 kpc reaches ~ 10^ K, or 
even higher. The H2 and HD molecule abundances are 
consequently suppressed via effective collisional dissoci- 
ation, leading to values of ~ 10~^ and /hd ^ 10~^^, 
respectively. The mean density of the gas in the central 
100 pc is nn ~ 0.01 cm~^ which is 4 orders of magnitude 
smaller than that in the BHN run, showing the strong 
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Fig. 13. — I mpac t of HMXB feedback. We compare the protogalactic gas properties in simulations BHN (left) and BHB (right) at z = 15, 
similar to Fig. |11| but now extending to 30 kpc to illustrate the HMXB feedback effect not only on the gas within the host halo, but also 
on the extended IGM. 



photoevaporation due to HMXB radiative feedback. 

In both simulations part of the IGM outside a 1 kpc 
radius is heated and ionized by an I-front which had been 
generated 6.5 Myr ago by Pop III stars formed ~ 5 kpc 
away from the center of the halo hosting the HMXB. 
The HMXB feedback effect on the IGM is most clearly 
evident in gas which has never experienced any stellar 
radiation feedback, as the latter could masquerade as 
the former. The distant IGM in the BHN simulation 
is in a cold and neutral state, while in the BHB run, 
the IGM unaffected by any stellar radiative feedback is 
in a warm and partially ionized phase because of the 
pervasive effect of the strong HMXB emission. These are 
favorable conditions for H2 and HD molecule formation, 
promoting their abundance by an order of magnitude 
over the BHN case. 



To more clearly distinguish the HMXB feedback on 
the first galaxy and on the distant IGM, we compare 
in Figure [1^ the volume- averaged temperature, electron 
fraction, and H2 fraction for simulations BHN and BHB 
inside and outside of a 15 kpc radius, beyond which our 
algorithm does not trigger the build-up of H II regions 
around newly-formed Pop HI stars. Hence, any varia- 
tions in gas properties in the far-zone are caused by the 
HMXB feedback alone. Here, we adopt a 15 kpc cutoff 
radius rather than 10 kpc, the region insi de w hich stars 
are endowed with an H II region (see Sec. 2.4), because 
an ionized region generated by a Pop III star formed near 
r ~ 10 kpc from the center of the emerging galaxy will ex- 
tend out to ~ 15 kpc. The corresponding X-ray heating 
increases the volume- averaged temperature within 15 kpc 
by only a factor of 3, while the temperature in the dis- 
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Fig. 14. Impact of HMXB feedback on IGM. We show the 
volume-averaged temperature, electron fraction, and H2 fraction 
vs. redshift in simulations BHN and BHB. We show these quanti- 
ties separately for radial distances smaller and larger than 15 kpc, 
beyond which the feedback from Pop HI stars is not included. 
Hence, the feedback effect from the HMXB on the distant IGM 
becomes clearly apparent. 

tant IGM is enhanced by two orders of magnitude. The 
difference in the near zone is relatively small because the 
gas within the 15 kpc radius is affected by both the feed- 
back from Pop III stars, and from the HMXB, indicating 
that the overall gas properties are dominated by the lo- 
cal stellar feedback. In the absence of any local stellar 
feedback, as is the case in the region beyond the 15 kpc 
radius, on the other hand, the IGM is substantially pre- 
heated, and partially ionized when HMXB feedback is 
present. These conditions facilitate molecule formation, 
enhancing abundances by an order of magnitude. Our 



results are in go od agreement with the study of Kuhlen 
& Madau (2005), given that the X-ray flux in our sim- 
ulation BHB is comparable to that for their power-law 
case. 

We assume that the HMXB source has been turned 
on for ~ 150 Myr, which is much longer than any real- 
istic lifetimes. For example, the main-sequence lifetime 
of a donor star of SOM^;^, close to the predicted typi- 
cal Pop HI mass (Clark et al.l 2Qllb[ Hosokawa et al. 



2QTT| [Stacy et al.H:^Ul2), would be ^ 20 Myr. Thus 



our work can be considered to represent the case where 
at least one HMXB exists continuously over 150 Myr, 
somewhere in the Lagrangian volume of the emerging 
protogalaxy. A preheated IGM would provide us with 
a promising observational window into the early uni- 
verse. Under these conditions, the IGM would be suf- 



ficiently hot to be decoupled from the CMB tempera- 
ture, TcMB = 2.73 K(l + z), but not hot enough to be 
ionized, leaving the IGM substantially neutral. Such a 
preheated, neutral IGM is likely detectable in the emis- 
sion of redshifted 21cm radiation, which is produced by 
the spin-flip transition between the singlet and triplet 
hyperfine- structure levels of neutral hydrogen. In prin- 
ciple, future 21 cm observations could provide constraints 
on the character of early feedback processes, and in par- 
ticular on the importa nce of a co ntribution from Pop HI 
miniquasars (see e.g. Furlanetto et al. 2006) . 

4.3. Comparison to SN Feedback 

In the previous section we discussed the situation es- 
tablished under radiative feedback from an accreting iso- 
lated BH, and an HMXB. Effectively, we have made the 
assumption that all Pop HI stars end their lives as a 
BH, without any SN explosion preceding their deaths. 
If a Pop HI star, however, has a mass in the range of 
140 — 260 M0, it would die as a PISN, exerting a strong 
mechanical feedback on the surrounding primordial gas, 
enriching it with the metals produced before and dur- 
ing the explosion. We are able to dir ectly compare our 
results to those of Greif et al. (2010), where PISN en- 
richment of the IGM, stellar radiative feedback, chemical 
mixing, and metal-line cooling were included, because all 
simulations start from the same initial conditions. The 
underlying DM structure was thus identical, providing 
us with an ideal laboratory to study the impact of varia- 
tions in the baryonic physics on the assembly process of 
the first galaxies. 

We find that the gas properties within the virial ra- 
dius of the host halo in our BHN simulation at z ~ 10 
are very similar to those in the Gre if et al.| (2010) PISN 
simulation, except that the central gas in tne latter sim- 
ulation is enriched to Z ^ IO^^Zq. It is, however, still 
unclear whether such enrichment level is sufficient to en- 
able the transition from Pop HI to Pop II star formation, 
because the critical metallicity threshold for this transi- 
tion is still a matter of debate. The additional heating 
from the central BH in our BHS run, on the other hand, 
significantly enhances molecule formation, possibly by up 
to two orders of magnitude, and the gas might be able to 
further fragment. Which feedback is dominant in shap- 
ing second- generation star formation in the first galax- 
ies? The answer will strongly depend on the IMF of the 
first generation of stars, and on the detailed physics of the 
transition from the predicted top-heavy IMF for Pop HI 
to the more normal IMF in already metal-enriched gas. 

5. SUMMARY AND CONCLUSIONS 

We have studied how the assembly of a primordial 
galaxy is affected by the radiative feedback from an 
accreting, isolated stellar-mass black hole and a high- 
mass X-ray binary (HMXB), which are two possible 
end products of Pop HI star formation. To accomplish 
this, we have carried out three cosmological simulations 
which self-consistently account for the radiation from 
individual Pop HI stars, and from a central black hole 
X-ray source. In particular, we have focused on the 
early growth of a Pop HI black hole remnant at high 
redshift s, and on the role of black hole feedback in 
shaping subsequent star formation in the emerging first 
galaxy at z ~ 10. Our main results and implications can 
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be summarized as follows. 

1. We have shown that locally, as opposed to any 
global effect on the large-scale IGM, the feedback from 
an isolated, accreting black hole is very efficient, lead- 
ing to a strong suppression of the early growth of the 
seed black hole. Without such feedback, the growth 
rate quickly reaches near-Eddington values. Our results 
imply the following important consequence. A widely 
adopted assumption for the formation of early super- 
massive black holes (SMBHs), whose presence is inferred 
from QSO observations at z ~ 6, is that their growth 
begins with a Pop III seed black hole. If this assumption 
were correct, quenching of black hole accretion would 
be required to av oid an overabundance of SMBHs at 
low re dshifts (e.g., |Bromm fc Loeb||2QQ3| |Trenti fc Sti-| 
2QQ7| ). Without such quenching, one^ould also 



avelli 



run mto mconsistencies with the low number density 
of quasars, (6 ± 2) x 10~^^Mpc~^ at z ~ 6, assu ming 
Ho = 65km s~^Mpc~\l],^, 0.35, and 0.65 ( |Fan 



et al.|2 QQ4l). We suggest that the radiative feedback from 



accreting black holes plays a key role in suppressing early 
black hole growth, thus constraining models for SMBH 
formation. 

However, in this study we have tracked only an iso- 
lated seed black hole, born in a halo of 5 x 10^ Mq that 
virialized at z = 28. Alternatively, a black hole seed 
could form significantly earlier, in a halo corresponding 
to rare, very high-<j peaks, thus providing a head start to 
the build-up process. A differe nt growth-boosting mech - 
anism has been suggested by |Volonteri & Rees ( 2005 ) , 
arguing that Pop III remnant black holes, hosted in ha- 
los with a virial temperature of Tvir ^ 10^ K at z ~ 25, 
corresponding to peaks of 4cr or above, may be able to 
experience stable supercritical accretion owing to effec- 
tive atomic hydrogen line cooling (see also Li|2011 ). It is 
still challenging to numerically test this scenario because 
of the prohibitive computational expense to capture such 
rare massive halos and follow their evolution. 

2. Despite the extremely small accretion rates onto the 
isolated black hole, the corresponding radiative feedback 
on the thermal and chemical evolution of the gas in the 
immediate surroundings, at scales < 100 pc, is signifi- 
cant. The black hole feedback establishes two distinct 
thermodynamic branches in the central region of the 
emerging galaxy: A hot phase with T ~ 10^ K and high 
ionization fraction, as well as a warm one with T ~ 10^ K 
and partially ionized conditions. In the latter phase, 
molecular hydrogen formation is favored, leading to an 
elevated fraction of fu^ > 10~^, two orders of magnitude 
higher than in the absence of the radiation from a black 
hole. This enhanced H2 fraction could promote second 
generation star formation within the emerging first galax- 
ies, possibly resulting in a burst-like mode. However, 
we do not include the effect of Lyman- Werner molecule- 
dissociating photons emitted by the soft (MOD) com- 
ponent of the accreting black hole. Our H2 and HD 
molecule abundances are therefore upper limits. In fu- 
ture work, we will explore the impact of Lyman- Werner 
radiation generated by black hole accretion. 

What would be the characteristics of the stars formed 
out of such metal-free gas with a significantly enhanced 
electron fraction, the so-called Pop HI. 2 stars (jO^Shea 



et al.||2008D . This issue is far from being settled. Exist- 
mg studies suggest that the elevated H2 and HD fractions 
will lead to rapid cooling, resulting in smaller mass stars, 
of the order of ^ 10 M©, compared to the canonical value 
of ~ 100 Mq fo r Pop II I.l stars formed in quasi-neutral 



minihalos (e.j 



O'Shea et al. 2005 Johnson & Bromm 



2QQ6l|Yoshida et"arpnD7b[). h owever, the suite of sim- 



ulations carried out by 'Clark et al. (2011a), who con- 
sidered turbulent fragmentation m primordial gas, has 
indicated, somewhat counter-intuitively, that the second 
generation. Pop HI. 2 stars were likely more massive than 
the Pop III.l stars, formed in the pristine minihalos that 
are unaffected by prior star formation. This can be un- 
derstood by the inability of the gas in the Pop HI. 2 case 
to reach the CMB temperature floor via HD cooling un- 
der realistic halo conditions. 

3. The feedback from an efficiently radiating HMXB is 
very strong locally, and moderately important globally. 
Locally, the effect on the surrounding primordial gas is 
to heat it to high temperatures of > lO^K, and to fully 
ionize it. The corresponding strong photo-evaporative 
outflow suppresses central gas densities, thus prevent- 
ing any subsequent star formation within the emerging 
galaxy. Our results imply that once a halo of ~ 10^ Mq 
harbors an HMXB, the ensuing strong radiative feedback 
will delay the condensation of gas in the atomic cooling 
halo, possibly leading to a decrease in the number of flrst 
galaxies at a given epoch. While we explore in this pa- 
per the impact of a continuously radiating HMXB, the 
strength of the effect will depend on the realistic abun- 
dance and lifetime of such HMXB sources at high red- 
shifts. The highly penetrating X-ray photons emitted 
by the HMXB are able to effectively heat and partially 
ionize the distant IGM as well, and as a result catalyze 
the formation of molecular hydrogen over large volumes. 
We have shown that the global star formation rate in the 
presence of HMXB feedback is, at least at early times, 
higher than in the runs without it. We assume that 
the broad protostellar mass function seen in recent high- 
resolution simulations suggesting that primordial stars 
have less extreme masses is likely to increase the pos- 
sibility of HMXBs. However, the mapping between the 
protostellar fragments, as represented by sink particles, 
to the flnal Pop III stars, is still incompletely understood. 
Improved simulations sh ould soon be able to further elu- 
cidate this key issue (see'Greif et al.| |2012D . 

To sum up, the radiative feedback from accreting 
black hole sources, either isolated or in a binary system, 
has a dramatic impact on the surrounding gas within 
the flrst galaxy. The alternative, where no black hole 
feedback is present, but where the gas is processed by 
a preceding PISN explosion, is less affected in the long 
run. However, it is likely that the enhanced cooling 
due to metals and dust will only become effective at 
densities that are higher than can be resolved in these 
simulations. One possible consequence of metal and dust 
cooling might be an enhanced star formation efflciency, 
potentially enabling the formation of a bona fide P op II 



stellar cluster inside the first galaxi es (Clark et al.pOOS^ 
Tsuribe fc Omukai||2008l [Dopcke et al. 2011):T^etter 



understand such Fop II star formation, we need to follow 
the subsequent collapse of the dense, metal-enriched 
gas with much higher numerical resolution in future 
simulations. 
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